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ABSTRACT 

The recent discovery of Jupiter-mass planets orbiting at a few AU from their 
stars compliments earlier detections of massive planets on very small orbits. The 
short period orbits strongly suggest that planet migration has occurred, with the 
likely mechanism being tidal interactions between the planets and the gas disks 
out of which they formed. The newly discovered long period planets, together 
with the gas giant planets in our solar system, show that migration is either 
absent or rapidly halted in at least some systems. We propose a mechanism 
for halting type-II migration at several AU in a gas disk. Photoevaporation of 
the disk by irradiation from the central star can produce a gap in the disk at 
a few AU, preventing planets outside the gap from migrating down to the star. 
This would result in an excess of systems with planets at or just outside the 
photoevaporation radius. 
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1. Introduction 

It is believed that gas giants do not generally form at small orbital distances from the 
central star (Boss 1995). Thus, a natural explanation for extra solar planets orbiting close 
to the central star is that these planets formed further away in the protoplanetary disk 
and migrated inward to where they are now observed. A variety of mechanisms have been 
proposed to explain planet migration: the interaction between a planet and a planetesimal 
disk (Murray et al. 1998), the gravitational interaction between two or more Jupiter mass 
planets (Rasio & Ford 1996), and the tidal gravitational interaction between the planet and 
the surrounding disk gas (Goldreich & Tremaine 1979, 1980). The last mechanism, focused 
on in this paper, is expected to be dominant at early times since the surrounding gaseous 
disk is required for the formation of planets. 

If the perturbation exerted on the disk by the planet is small, the disk structure is not 
greatly altered, and the planet moves inward relative to the surrounding gas (Ward 1997). 
This type of migration is referred to as 'type-I'. However, if the planet is large it may open 
a gap in the disk (Goldreich & Tremaine 1980). The planet is locked to the disk and moves 
either inward or outward in lock-step with the gaseous disk. This slower migration is referred 
to as 'type-II'. 

We propose a mechanism for halting type-II migration: photoevaporation driven by 
radiation from the central star. The planet's final location is consistent with the Solar Sys- 
tem and the growing class of extra solar planets with nearly circular orbits outside of a few 
AU (see Tinney et al. 2002, figure 4). Photoevaporation by the central star was proposed 
by Shu et al. (1993) and Hollenbach et al. (1994) as a way to remove a gas disk. Hollen- 
bach et al. (2000) generalized the discussion, describing the variety of possible disk removal 
mechanisms: accretion, planet formation, stellar encounters, stellar winds or disk winds, and 
photoevaporation by ultraviolet photons from the central source or massive external stars. 
Hollenbach et al. (2000) concluded that the dominant mechanisms for a wide range of disk 
sizes are viscous accretion and photoevaporation, operating in concert within the disk. In 
this paper, we consider photoevaporation by the central source and viscous accretion. 

2. Model 

The model for disk removal used here is developed in a paper by Matsuyama, John- 
stone, & Hartmann (2003) and is similar to that used by Clarke, Gendrin, & Sotomayor 
(2001). Additionally, we assume a planet with a large mass, which opens a narrow gap in 
the disk, and assume that planet migration proceeds in lock-step with the disk evolution 
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(type-II migration). The gas disk orbiting the central star with the local Keplerian circular 
velocity, Vk, is axisymmetric and geometrically thin. Considering angular momentum and 
mass conservation of a disk annulus at a radius, R, with kinematic viscosity, is, the disk 
surface density evolution can be described by (Pringle 1981) 



<9£ _ 3 d 
~dt ~ RdR 



Rl/2 M 



where £ is the disk surface density and t is the disk evolutionary time. We adopt the 
standard a-parameterization of Shakura & Sunyaev (1973) and write v = ac s H, where c s is 
the sound speed at the disk mid-plane and H is the disk thickness. Hartmann et al. (1998) 
estimate a ~ 1(T 2 - 1(T 3 and a disk mass ~ 0.01 - O.2M for T Tauri stars (TTSs). For 
the modeled disk temperature distribution, oc R~ 1 / 2 (D'Alessio et al. 1998), the viscosity 
takes the simple form, v(R) oc R. Given a solution for equation 1, we find the drift velocity, 

and describe the evolution of the disk stream lines. 

The EUV (A < 912A) photons from the central star and the accretion shock (we will 
refer to these photons as the EUV photons from the central source) are capable of ionizing 
hydrogen and evaporating material from the disk surface. Photoevaporation forms an ionized 
atmosphere above the thin viscous disk. Hollenbach et al. (1994) found analytic solutions for 
the photoevaporation mass loss rate by EUV photons from the central source. These photons 
are attenuated by recombined hydrogen atoms and scattered in the ionized atmosphere, 
providing a source of diffuse EUV photons that dominates the flux onto the disk at radii 
much larger than the size of the central star. Disk material is gravitationally bound to the 
central star inside the gravitational radius, R g = GM*/3cf, and it flows out of the disk 
surface at the sound speed of the heated material, c' s , beyond this radius. Given the number 
density of ionized hydrogen at the base of the ionized layer, n(R), we can calculate the 
evaporation rate: 

£ (r) = 1 ty™ p n(R)d a , if R > Rg; ( o\ 

ph{ ' \0, otherwise, 1 ' 

where \i = 0.68 is the mean molecular weight for the ionized material. 

Assuming ionization equilibrium Hollenbach et al. (1994) found the number density at 
the base of the ionized layer, for R > R g : 

»'^^ llS (it) 1,! (a)(ii)" ! -" !: (4) 
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where <fi is the total ionizing flux from the central source. We assume that half of the accretion 
luminosity is radiated as hot continuum from accretion shocks and that the EUV flux can 
be characterized as blackbody emission at T as = 10 4 K (Gullbring et al. 2000; Johns-Krull 
et al. 2000; Kenyon et al. 1989). There are currently no strong observational constraints on 
the ionizing flux from the central star, 0*. The EUV flux for the sun is ~ 2.6 x 10 37 s _1 
(see Wilhelm et al. 1998, figure 10b). Since TTSs are chromospherically active (D'Alessio 
et al. 1998), their EUV fluxes must be higher by few orders of magnitude. We estimate 
the ionizing flux from the FUV spectra of TTSs (Valenti et al. 2000, figure 3) and the 
corresponding distances to the star (Johns-Krull et al. 2000, table 1). The typical values are 
~ 10 39 s _1 for classical T Tauri stars. Extrapolating the FUV flux underestimates the EUV 
flux; thus the EUV flux might be as high as lO^s -1 . 

Combining photoevaporation with viscous accretion is done numerically. At each time 
step, photoevaporation induced mass loss and viscous diffusion induced disk evolution are 
solved. For details on the disk dispersal calculation refer to Matsuyama et al. (2003). The 
simulation is stopped when the disk mass is IMj. 

3. Discussion 

At a given disk radius, a gap forms when the mass transported by viscous accretion is 
equal to the mass removed by photoevaporation. This occurs when the photoevaporation 
time scale, E/E p ^, is equal to the viscous diffusion time scale, R 2 /3u. Therefore, we can 
estimate the surface density for gap formation at the gravitational radius, 

^M2gc m -)( I ii 5 )- 1 ( I ^) 1/2 (i)", (5) 

Figure 1 shows snapshots of the disk surface density distribution for a disk with initial 
mass, M(0) = 0.03M© and a = 10~ 3 . A gap starts forming at the gravitational radius 
when the surface density reaches S 9ap ~ 2 gem" 2 during the last stages (t gap ~4x 10 7 yr) 
of the disk evolution and the disk is divided into two annuli. The subsequent evolution 
is dominated by two competing effects; viscous diffusion attempts to spread both annuli 
and remove the gap structure while photoevaporation removes material predominantly at 
the gravitational radius, reopening the gap. The material in the inner annulus is quickly 
removed both at the inner edge, by accretion onto the central star, and at the outer edge, by 
the combination of viscous spreading and photoevaporation. In contrast, the outer annulus 
loses material primarily from its inner edge. The inner edge of the outer annulus recedes 
from the star as photoevaporation removes disk material, thereby reducing the mass removal 
rate in accordance with equations (3) and (4). 
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We model a planet which has opened a narrow gap in the disk. The migrating planet 
will open a gap if its mass is greater than the critical mass, M c ~ M* (81nacj./8vl), (see, e.g., 
Lin & Papaloizou 1993). The planet is initially tidally locked to the disk and moves with 
the gas, migrating outward or inward depending on its initial location. This migration may 
stop for two reasons. First, the combination of photoevaporation from the central source 
and viscous spreading removes disk mass, reducing the gravitational interaction between 
the planet and the disk and allowing the planet to decouple its orbit from that of the 
gas. Second, the planet will stop if it reaches the gap created by photoevaporation, as the 
tidal torques are unable to influence the planet across such a large opening. In order to 
describe planet migration, we show in figure 2 the disk stream lines which start between 4 
and 18AU, at t — 10 4 yr, where there is enough mass to form a Jupiter mass planet in an 
annulus of width 4AU. After the planet opens a gap it migrates along the disk stream lines; 
therefore, we can follow planet migration in figure 2 for different planet gap opening locations 
and times. The disk mass outside the gravitational radius is reduced to Jupiter's mass in 
~ 10 8 yr (see figure 3); therefore, a migrating planet with Jupiter's mass stops migrating at 
any disk location soon after this time. However, a migrating planet may reach the inner disk 
boundary (R = 1CT 2 AU) or the photoevaporation gap before the disk mass is reduced to 
Jupiter's mass. 

Figure 2 provides information on which initial locations, and times, allow a planet 
tidally locked to the disk to survive. It is difficult to time the gaseous disk removal in such 
a manner as to maintain planets at large distances from the central star. A planet which 
opens a narrow gap at t — 10 4 yr only survives if its initial disk location is greater than the 
critical radius ~ 14 AU, because the planet reaches the outer edge of the photoevaporation 
gap after the gap has formed. Since the disk expands rapidly at early times, this critical 
radius increases with later planet gap opening timescales; i.e. the later the planet forms, the 
further away from the central star it has to be in order to survive. However, this behavior 
changes after the photoevaporation induced gap forms. If the planet forms at a time similar 
to when the photoevaporation gap forms, t gap ~4x 10 7 , it survives as long as it originates 
outside the gravitational radius. 



4. Conclusions 

In this letter, we have studied different scenarios in which a planet undergoing type- 
II migration survives due to the formation of a gap by photoevaporation from the central 
source. If the planet gap opening occurs close to when the photoevaporation gap forms (at 
i~4x 10 7 yr for our fiducial model), the final semi- major axis of the planet is near the 
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gravitational radius (~ 2.4AU) for a wide range of initial semi-major axes. This scenario 
requires long time scales for planet formation (10 6 — 10 8 yr), similar to those predicted by 
the core accretion model (Pollack et al. 1996; Bodenheimer et al. 2000). 

We can also predict the final semi- major axis when considering rapid planet formation 
by disk instability (Boss 2000). Our disk initially spreads outward under the influence of 
viscosity, allowing planets that form early at small radii to move out beyond R g before 
the photoevaporation gap opens, and enhancing their survival rate. The rate and amount 
of viscous spreading in real disks depends on the initial surface density profile and any 
continuing accretion onto the disk. If the planet induced gap opening occurs at early stages of 
the disk evolution (for example, at t ~ 10 4 yr), all the planets initially located inside ~ 14AU 
will reach the disk inner boundary at the central star. The only planets which survive are 
those which form at a narrow range of semi-major axes (14 < R < 18AU). However, viscous 
diffusion spreads this narrow range during planet migration, and the final semi-major axis 
distribution is 3 < R < 22AU. This final distribution of planets has similarities to our 
own Solar System and it suggests that the gas giants could have formed at ~ 14 — 18 AU 
and migrated to their current locations. Furthermore, the initial planet locations are in 
agreement with recent simulations of planet formation by disk fragmentation (Mayer et al. 
2002). 

The minimum surface density for gap opening by photoevaporation, S ffap , and the disk 
size, Rd, determine the disk mass outside the gravitational radius. In our fiducial model, the 
mass outside the gravitational radius at the time of gap opening is 27rS gapRgi^Rd — ^g) ~ 
3Mj; therefore it is possible to consider migration of a Jupiter mass planet. However, for 
lower ionizing fluxes (i.e. <fi± < 10 40 s~ 1 ), the minimum surface density for gap opening 
by photoevaporation becomes smaller (see equation 5) and the halting of planet migration 
becomes marginal. For example, for 0* = 10 38 s _1 , the minimum surface density for gap 
opening and the disk mass outside the gravitational radius are smaller by one order of 
magnitude (S 9ap ~ 0.2gcm~ 2 and M out = O.SMj). 

We note that the gravitational radius, R g , roughly coincides with the location of the 
asteroid belt in the solar system. It is well known that the surface density of heavy elements 
reaches a severe local minimum in the asteroid belt; this is usually attributed to long term 
gravitational effects from Jupiter, although it is not clear that the amount of depletion that 
results is sufficient (Liou & Malhotra 1997; Holman & Murray 1996). Simple estimates show 
that the planetesimal formation time scale is much shorter than the gap formation time, so 
it is not clear that the formation of a photoevaporation gap would significantly reduce the 
surface density of planetesimals, but the coincidence is intriguing. 

The photoevaporation gap opening timescale (~ 4 x 10 7 yr) is long compared to the 
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typically observed gaseous disk lifetimes (10 6 — 10 7 yr). However, the disk lifetime for stars 
outside stellar clusters is considerably longer than the disk lifetime for stars in the hostile 
stellar clusters (Matsuyama et al. 2003). In addition, there is evidence that some stars retain 
disks between 10 7 — 10 8 yr (see Hillenbrand 2002, figure 2). 

It is also possible to halt both type-I and type-II migration in stellar clusters by removing 
the disk. This is a likely situation for the disks surrounded by ionization fronts in the 
Trapezium Cluster (Johnstone et al. 1998; Storzer & Hollenbach 1999; Bally et al. 2000). 
Storzer & Hollenbach (1999), Hollenbach et al. (2000) and Matsuyama et al. (2003) show 
that the disk lifetime due to photoevaporation by external stars is in the range 10 5 — 10 6 yr 
in the neighborhood of massive O stars, and between 10 6 — 10 7 yr at reasonable distances 
(~ 0.03 — 0.3pc) from external O stars. These timescales are in agreement with observational 
estimates for the lifetime of protoplanetary disks in stellar clusters (Haisch et al. 2001). Since 
the outer disk is removed first, there is no photoevaporation gap formation. This scenario is 
possible when considering the rapid planet formation by disk instability; on the other hand, 
it is uncertain whether gas giant planets form by core accretion in these short lived disks. 
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Fig. 1. — Snapshots of the surface density for a fiducial model under the influence of viscous 
diffusion and photoevaporation from the central source. The model corresponds to a — 10~ 3 , 
and an initial disk mass, M(0) = 0.03M Q . The short-dashed line indicates the location of the 
gravitational radius and the long-dashed line corresponds to the minimum surface density 
for gap formation by photoevaporation. The solid curves represent t — 0, 10 6 , 10 7 , 3.6 x 
10 7 ,4.2 x 10 7 ,4.7 x 10 7 ,5.2 x 10 7 ,5.7 x 10 7 ,6.2 x 10 7 ,6.7 x 10 7 , 7.2 x 10 7 , 7.7 x 10 7 ,8.2 x 
10 7 , 8.7 x 10 7 , and9.2 x 10 7 yr. The gap structure starts forming at t gap ~4x 10 7 yr, when the 
disk mass is ~ 3 x Mj and the surface density at the gravitational radius is S 9ap ~ 2g cm~ 2 
(long-dashed line). 
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Fig. 2. — Disk stream lines for the fiducial model with M(0) = O.O3M and a = 10~ 3 . The 
initial spacing (at t = 1 4 yr) between the streamlines is 1AU, covering the disk from 4AU to 
18AU. The short-dashed line indicates the location of the gravitational radius, R g = 2.4AU, 
and the dotted line represents the time when the photoevaporation gap starts forming, t gap ~ 
4 x 10 7 yr. Planets forming at R ~ 5AU and tidally locking to the disk early (t = 10 4 yr) 
(open triangles) or late (solid triangles) fall onto the central star at t ~ 2.5 x 10 6 yr. In 
contrast, a planet starting with initial semi-major axis R ~ 15AU at t — 10 4 yr (crosses) 
stops migrating when it reaches the photoevaporation gap. 
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Fig. 3. — Disk mass outside the gravitational radius, M out , and disk mass inside the grav- 
itational radius, M in , as a function of disk lifetime for the fiducial model. The dotted line 
indicates the time, t gap ~4 x 10 7 yr, when the gap structure starts forming. 



